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ABSTRACT 

The  Naval  Oceanographic  Office  developed  the  Integrated  Carrier 
Acoustic  Prediction  System  (ICAPS)  in  order  to  perform  on-scene 
acoustic  forecasts  for  the  Fleet.   To  evaluate  the  system,  bathy- 
thermographic  information  was  obtained  for  points  in  the  North 
Atlantic  during  the  naval  exercise  "SEACONEX."   Fleet  Numerical 
Weather  Central  (FNWC)  provided  climatological  and  synoptic  temx^era- 
ture  profiles  along  with  related  acoustical  forecasts  for  the  points 
under  study.   Similarly,  ICAPS  supplied  a  combined  climatological 
and  synoptic  temperature  profile  and  related  acoustic  forecasts  for 
each  of  the  points.   Comparisons  of  the  two  systems  were  made 
against  an  actual  mean  temperature  profile  and  the  acoustic  fore- 
cast that  it  generated.   It  appears  that  the  climatological  acoustic 
forecasts  were  not  reliable  tactically.   Yet,  the  synoptic  forecasts 
by  FNWC  did  yield  reliable  low  frequency  acoustic  forecasts  within 
the  limits  of  computer  roundoff  error.   ICAPS  capability  for  only 
a  single  temperature  profile  input  proved  to  yield  reliable  acoustic 
forecasts,  and  should  be  used  as  the  real-time  acoustic  forecast 
system  for  the  Fleet. 
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I.   INTRODUCTION 

A.   OBJECTIVES 

The  primary  objective  of  this  study  is  to  compare  the  ocean 
thermal  structure  analyses  and  acoustical  products  of  Fleet 
Numerical  Weather  Central  (FNWC)  with  the  Integrated  Acoustic  Pre- 
diction System  (ICAPS)  thermal  structure  inputs  and  acoustic  fore- 
casts.  This  comparison  is  intended  to  evaluate  the  overall 
credibility  of  both  systems.   The  data  used  in  this  study  were 
collected  aboard  the  USS  GUAM  (LPH-9)  during  the  naval  exercise 
"SEACONEX"  held  in  the  North  Atlantic  Ocean  in  June  1973. 

In  order  to  accomplish  this  objective  (basic  plan  including 
discussion  of  FNWC  models,  etc.)  this  paper  presents  a  brief  dis- 
cussion of  the  merits  of  the  climatology  used  by  FNWC  and  the  new 
interim  climatology  developed  by  the  Naval  Oceanographic  Office 
for  the  ICAPS  system.   Additionally  the  FNWC  synoptic  temperature 
profile  and  its  related  acoustic  field  are  compared  to  those 
generated  by  using  average  mean  bathythermograph  (BT)  profiles 
obtained  at  sea.   The  value  of  a  single  BT  in  the  real  ocean  is 
discussed.   In  the  case  of  a  single  BT  profile  obtained  aboard 
ship,  a  theory  concerning  the  errors  associated  with  using  such 
a  profile  is  advanced. 

A  second  objective  of  this  study  is  to  appraise  the  operational 
techniques  required  for  ICAPS.   By  a  detailed  discussion  of  the 
strengths  and  weaknesses  of  ICAPS  temporary  profile  inputs  combined 


with  suggestions  on  how  to  capitalize  on  the  strengths  and  compen- 
sate for  the  weaknesses,  a  pattern  of  operation  will  be  projected. 
Planning  and  operating  procedures  for  the  use  of  ICAPS  are  proposed. 
Discussions  are  aimed  at  the  procurement  and  processing  of  oceano- 
graphic  data  rather  than  the  tactical  utilization  of  such  informa- 
tion. The  goal  of  this  portion  of  the  study  is  to  improve  the 
effectiveness  of  the  shipboard  acoustical  forecaster. 

B.   ACOUSTIC  FORECASTING  MODELS 

The  Fleet  Numerical  Weather  Central  (FNWC)  in  Monterey,  Cali- 
fornia, currently  produces  and  distributes  to  the  fleet  two  primary 
acoustic  forecasting  products,  the  Ship,  Helicopter  Acoustic  Range 
Prediction  System  (SHARPS  II)  and  the  Acoustic  Sensor  Range  Predic- 
tion (ASRAP  III)  (Navweaservcom,  1971a, b)  .. 

1 .   Ship,  Helicopter  Acoustic  Range  Prediction  System  (SHARPS  II) 
The  Ship,  Helicopter  Acoustic  Prediction  System  (SHARPS)  is 
a  computerized  sonar  range  prediction  scheme.   It  is  designed  to 
take  environmental  analyses  and  convert  them  into  acoustic  range 
predictions.   Inputs-  to  the  program  include  average  equipment  para- 
meters, platform  speed,  sonar  operating  mode,  sensor  and  target 
depth.   Range  predictions  are  based  upon  a  single-ping,  50  percent 
probability  of  target  detection  by  an  unalerted  sonar  operator. 

The  key  input  to  the  SHARPS  II  program  is  the  vertical  tem- 
perature profile o   This,  along  with  the  climatological  salinity 
profile,  is  transformed  by  the  use  of  Wilson's  (1960)  equation 
for  sound  speed  into  a  sound  velocity  profile  (SVP) .   In  turn,  the 
SVP  is  used  to  generate  a  raytrace  pattern  considering  direct  path, 


purely  refracted,  and  the  surface  and  bottom  reflected  propagation 
paths.   The  change  in  sound  intensity  with  range  is  computed  con- 
sidering losses  due  to  spreading,  scattering,  reverberation,  absorp- 
tion, and  surface  and  bottom  reflections.   Average  performance 
characteristics  are  used  for  each  sonar  type,  and  an  estimate  of 
the  sound  intensity  as  a  function  of  range  is  achieved.   Considering 
that  for  active  sonars  a  two-way  travel  distance  (distance  from 
source  to  target  and  return) is  needed;  the  amount  of  propagation 
loss  is  double  that  of  the  one-way  travel. 

SHARPS  II  claims  increased  accuracy  over  previous  models. 
This  accuracy  is  attributed  to  (1)  use  of  the  modified  Atlantic 
Meteorological  and  Oceanographic  Study  (AMOS)  equations  for  surface 
channel,  direct  path,  and  sub-surface  duct  computations;  (2)  inclu- 
sion of  the  phys.ics  of  the  operational  raytrace  and  propagation 
loss  model  techniques  for  bottom  bounce  and  convergence  zone  compu- 
tations; (3)  inclusion  of  the  basic  reverberation  technique  used  in 
the  Navy  Interim  Surface  Ship  Sonar  Prediction  Model  (NISSM) 
developed  by  Mobile  Sonar  Technology  (MOST)  Committee;  and  (4)  in- 
creased input  data  to  the  model  including  a  description  of  the  SVP 
from  the  surface  to  the  bottom,  a  detailed  sonar  system  description 
for  all  equipments,  and  area  scattering  strengths  for  reverberation 
computations. 

Each  prediction  is  made  for  a  predetermined  acoustical  area 
in  which  several  assumptions  are  made:  (1)  water  depth  is  constant; 
(2)  SVP  is  horizontally  constant;  and  (3)  bottom  composition  and 
roughness  are  constant.   Predictions  are  made  on  the  basis  of  12 
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hourly  environmental  temperature  analyses  by  FNWC.   The  reliability 
and  variability  of  that  temj^erature  input  directly  affects  the 
reliability  and  variability  of  the  acoustic  forecast. 
2.   Acoustic  Sensor  Range  Prediction  (ASRAP  III) 

The  ASRAP  III  program  is  directed  toward  the  requirements 
of  the  airborne  acoustic  sensors,  primarily  the  passive  sonobuoys. 
Normally  two  distinct  forecasts  are  available:  a  forecast  of  passive 
propagation  loss  at  discrete  frequencies,  and  a  forecast  for  the 
active  airborne  sensors,  primarily  the  active  sonobuoys. 

ASRAP  III  uses  a  sophisticated  propagation  loss  model  which 
employs  an  analytical  solution  to  the  wave  equation  to  determine 
the  acoustic  intensity  for  various  ray  paths.   The  ocean  is  divided 
into  acoustically  homogeneous  areas  in  which  factors  that  affect 
sound  propagation  remain  constant  in  the  horizontal.   In  general, 
the  intensity  level  at  a  receiver  is  a  combination  of  sound  arriv- 
ing from  the  same  four  paths  mentioned  in  the  SHARPS  discussion. 
Losses  of  energy  by  reflection  from  a  boundary,  spreading,  and 
scattering  are  calculated  for  propagation  within  the  surface  duct 
and  below  it.   Of  course,  with  a  passive  system  only  one-way  trans- 
mission from  the  source  to  receiver  must  be  considered. 

The  SVP  used  by  the  ASRAP  program  is  derived  in  the  same 
manner  as  that  used  by  SHARPS.   Similarly,  the  reliability  of  the 
ASRAP  acoustical  forecast  is  a  function  of  the  reliability  of  the 
temperature  profile. 
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C.  INTEGRATED  CARRIER  ACOUSTIC  PREDICTION  SYSTEM  (ICAPS) 

The  Naval  Oceanographic  Office  (NAVOCEANO)  and  the  Naval  Under- 
sea Center,  San  Diego  (NUC)  /  have  developed  an  acoustic  prediction 
system  capable  of  providing  on-scene  forecasts.   A  mobile  unit 
capable  of  rapid  forecasts,  yet  small  enough  to  be  used  aboard  de- 
stroyers, submarines  and  VP  aircraft,  this  system  has  added  a  new 
capability  to  acoustic  prediction.   The  system,  now  called  the  In- 
tegrated Carrier  Acoustic  Prediction  System  (ICAPS),  is  composed 
of  a  Data  General  Corporation  Nova  800  computer,  a  Pertec  tape 
recorder,  a  Tektronix  4002A  graphic  computer  terminal  with  CRT, 
and  a  Tektronix  hard  copy  unit  (Consultec,  1972) . 

In  order  to  ensure  that  the  ICAPS  forecast  is  acoustically  valid, 
the  system  has  been  programmed  with  the  SHARPS  II  and  ASRAP  III 
models.   This  gives  the  unit  the  advantages  of  the  sophisticated 
acoustic  models  used  by  FNWC.   A  significant  difference  does  exist 
between  this  on-scene  approach  and  that  used  by  FNWC.   FNWC  uses 
available  BT  reports  from  the  preceding  72  hours  and  depends 
heavily  upon  climatology,  even  in  the  near  surface  region.   ICAPS 
uses  recent  BT  profiles,  or  composites  of  them,  in  the  near  surface 
region.   NAVOCEANO  has  developed  a  preliminary  climatology  for 
ICAPS  which  provides  the  thermal  structure  at  depths  greater  than 
range  of  the  BT. 

D.  TEMPERATURE  PROFILE  DERIVATION 

1.   Fleet  Numerical  Weather  Central 


FNWC  computes  the  thermal  structure  in  the  upper  1200  feet 
in  the  Northern  Hemisphere  every  12  hours.   This  analysis  is 
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composed  of  three  parts:  a  sea  surface  temperature  analysis,  a 
mixed  layer  depth  analysis,  and  a  sub-surface  temperature  analysis. 

a.  Sea  Surface  Temperature  Analysis 

The  present  scheme  used  by  FNWC  for  sea  surface  tempera- 
ture (SST)  analysis  is  the  Fields  by  Information  Blending  (FIB) 
technique  (Holl  et  al,  1971).   FIB  performs  numerical  analysis  of 
temperature  fields  by  the  assimilation  and  the  blending  of  informa- 
tion from  observations  and  climatological  averages  of  sea  surface 
temperatures.   It  combines  information  from  three  stages  in  time — 
current,  near  past  and  historical  temperatures.   Details  of  the 
analysis  procedure  are  given  in  Appendix  A. 

b.  Mixed  Layer  Depth  Analysis 

The  mixed  layer  depth  (MLD)  is  computed  by  the  POTMLD 
program,  which  takes  into  account  both  climatological  and  synoptic 
information  (Hesse,  1973) „   Details  are  given  in  Appendix  B. 

c.  Sub-surface  Temperature  Analysis 

The  sub-surface  analysis  model  as  described  in  the  new 
FNWC  Computer  Products  Manual  (Hesse,  1973)  is  constructed  on  the 
premise  that  the  number  of  BT's  available  in  any  one  day  is  small, 
much  too  small  to  make  a  truly  synoptic  analysis  possible.  Thus, 
it  seeks  to  make  use  of  the  BT  information  that  is  available,  but 
with  a  strong  dependence  upon  climatology  and  the  most  recent  SST 
analysis.   Details  are  given  in  Appendix  C. 

d.  FNWC  Climatology 

Stevens  (FNWC,  personal  communication)  indicated  that 
FNWC  developed  its  Atlantic  climatology  based  upon  National  Oceanic 
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Data  Center  (NODC)  station  data,  various  atlases  and  from  500,000 
-  1,000,000  mechanical  BT  observations.   This  data  was  weighted  so 
that  the  effects  of  a  small  amount  of  data  from  any  time  period 
would  not  be  treated  with  the  same  weight  as  a  larger  amount  of 
data  for  another  time  period  or  area.   This  step  was  necessary  to 
ensure  that  the  resulting  climatology  profile  would  be  representa- 
tive of  the  entire  time  period  to  which  it  pertained,  and  would  not 
be  dominated  by  any  limited  portion  of  the  data. 
2.   ICAPS  Temperature  Profile 

Fox  (NAVOCEANO,  personal  communication)  stated  that  either 
a  single  BT  observation  or  a  hand-prepared  composite  may  be  used 
as  an  input.   This  is  merged  at  1200  feet  by  the  computer  with  the 
ICAPS  climatology.   There  are  presently  no  capabilities  for  storage 
of  BT  reports  for  the  previous  72  hours  nor  is  there  a  function  for 
the  preparation  of  composite  BT  profiles  in  time  or  space.   Also, 
no  gross  error  check  (GEC)  is  performed  by  the  computer. 

a.   ICAPS  Climatology 

The  Atlantic  climatology  used  by  ICAPS  was  developed 
by  NAVOCEANO  as  an  interim  climatology  until  a  more  sophisticated 
model  could  be  developed.   Temperature  data  were  taken  from  NODC 
files  and  the  NAVOCEANO  oceanographic  atlas.   All  available  infor- 
mation within  each  five  degree  square  of  ocean  area  was  averaged 
at  standard  levels  to  obtain  the  climatological  profile.   An 
average  MLD  for  the  area  was  also  included  (Fox,  personal  communication) 
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II.   COMPARISON  OF  TEMPERATURE  AND  ACOUSTIC  PROFILES 

A.   BASIS  FOR  COMPARISONS 

The  purpose  of  this  portion  of  the  paper  is  to  compare  the  tem- 
perature profiles  used  by  FNWC  with  those  used  by  the  ICAPS  program 

r 

for  the  generation  of  their  respective  acoustical  forecasts.   As 
discussed  in  Chapter  I,  the  temperature  input  to  ICAPS  normally  con- 
sists of  the  latest  available  BT  to  a  depth  of  1200  feet.   Below 
1200  feet  the  BT  trace  is  merged  by  the  computer  with  the  ICAPS 
climatology  (Consultec,  1972) .   For  the  upper  1200  feet  FNWC  uses  the 
available  SST  and  BT  reports,  blending  them  in  the  horizontal  with 
climatology  and  the  analysis  for  the  previous  12  hours.   The  analyzed 
temperature  profile  is  then  merged  with  the  FNWC  climatology  below 
1200  feet  (Hesse,  1973).   When  BT  reports  are  not  available,  both 
FNWC  and  the  ICAPS  program  can  generate  acoustical  forecasts  using 
only  their  respective  climatologies. 

The  data  required  to  compare  the  FNWC  and  ICAPS  temperature  pro- 
files includes  on-scene  BT  reports,  the  FNWC  analyzed  thermal  struc- 
ture, and  the  FNWC  and  ICAPS  climatologies.   The  at-sea  data  were 
collected  on  board  USS  GUAM  (LPH-9)  in  the  North  Atlantic  Ocean 
during  June  1973.   These  profiles  were  obtained  by  the  ships  and 
helicopters  participating  in  exercise  "SEACONEX."   The  BT  profiles 
were  simultaneously  used  as  inputs  to  the  ICAPS  program  and  sent  to 
FNWC  as  inputs  to  their  thermal  structure  analysis  program.   The 
data  collection  points  are  indicated  in  Figure  1.   It  should  be  noted 
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that  due  to  the  nature  of  the  exercise,  considerably  more  data 
were  collected  at  points  B,  C,  and  D  than  at  the  other  points. 
While  on  board  USS  GUAM  the  ICAPS  climatology  profiles  were  also 
obtained.   Additionally,  during  June  the  FNWC  thermal  structure 
analysis  and  climatological  profiles  were  collected  and  held  at 
FNWC.   All  acoustic  forecasts  were  based  upon  the  same  salinity 
profile;  therefore,  that  parameter  does  not  enter  this  study  as 
a  variable. 

1.   Limitations  on  Data  Used  in  This  Study 

Admittedly,  because  of  the  severe  limitations  of  the  study 
with  regards  to  time  (limited  to  June)  and  space  (given  points  only) , 
conclusions  reached  must  be  treated  as  applicable  only  as  first 
approximations.   Another  constraint  which  was  present  during  this 
exercise  is  that  a  transiting  exercise  limits  the  number  of  BT's 
that  can  be  taken  in  an  area.   Thus  any  statistical  products  could 
be  in  error  on  the  same  order  of  magnitude  as  the  expected  varia- 
bility in  that  area.   Were  all  samples  to  be  taken  at  a  time  and 
place  experiencing  a  maximum  in  deviation  and  in  the  same  sense, 
then  a  biased  statistical  base  would  be  developed.   Thus  only 
definite  trends  in  this  study  will  be  regarded  as  indicative.   Extra- 
polations to  other  areas  or  time  periods  may  not  be  valid. 

B.   CLIMATOLOGY  COMPARISON 
1.   Temperature  Profile 

As  mentioned  in  Appendix  A,  at  FNWC  temperature  data  are 
analyzed  and  merged  at  designated  depths.   Since  most  BT's  obtained 
today  are  valid  at  least  to  a  depth  of  1200  feet,  all  comparisons 
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are  limited  to  these  levels0   The  basis  for  comparison  is  the  simple 
average  of  all  BT's  taken  at  a  position  during  a  specific  time  period; 
this  is  termed  the  actual  mean  temperature  profile.   Thus,  this  com- 
parison is  restricted  to  comparing  FNWC  climatologies  and  ICAPS  cli- 
matology for  the  upper  1200  feet.   As  differences  between  the 
climatologies  and  the  actual  mean  will  be  studied,  the  appreciable 
differences  will  be  assumed  to  occur  only  in  this  surface  area. 

Since  this  study  deals  with  deviations  from  a  mean  tempera- 
ture profile,  it  is  advantageous  to  separate  the  geographical  areas 
from  which  data  were  obtained  into  areas  where  the  vertical  tempera- 
ture profiles  exhibit  relatively  small  or  large  standard  deviations, 
respectively „   This  will  permit  an  investigation  of  any  relationship 
between  expected  standard  deviation  for  an  area  and  the  thermal  and 
acoustic  reliability  of  forecasts  for  that  area. 

To  properly  separate  the  positions  into  regions  of  high  or 
low  standard  deviation,  tables  of  expected  standard  deviation  con- 
structed by  FNWC  (Lewit,  1972)  are  used.   These  tables  cover  two 
degree  squares  over  all  ocean  areas.   Table  I  lists  the  values  for 
the  positions  used  in  this  study. 

Normalization  of  the  standard  deviation  for  each  depth  is 
accomplished  by  dividing  the  standard  deviation  at  that  depth  by 
the  maximum  standard  deviation  of  the  column.   These  results  are 
plotted  in  Figure  2  with  those  areas  judged  as  having  large  standard 
deviations  to  the  right  and  those  with  small  standard  deviations  to 
the  left.   This  judgment  is  based  on  the  requirement  that  at  800  feet 
the  curves  should  have  a  standard  deviation  less  than  0.3  and  at 
deeper  depths,  they  remain  below  0.3.   With  this  criteria,  artifically 
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Table  1. 


Position 

Depth 


Lewit's  (1972)  standard  deviation  (SD)  tables  for  positions 

of  int 

tions 

given 


SD 


tterest.   Included  under  total  is  the  number  of  observa- 
that  were  used  to  calculate  the  standard  deviation, 
in  degrees  centigrade. 

A  B  C  D  E 

Total  SD    Total  SD   Total  SD   Total  SD   Total 


,   00 

1.5 

56 

1.0 

57 

1.0 

57 

1.5 

10 

1.0 

4 

100 

l.l 

40 

1.0 

35 

1.0 

35 

0.9 

10 

0.8 

4 

200 

1.3 

39 

0.9 

35 

0.9 

35 

0.5 

10 

0.8 

2 

300 

1.1 

37 

0.8 

35 

0.8 

35 

0.6 

9 

0.9 

2 

400 

0.9 

32 

0.7 

34 

0.7 

34 

0.5 

9 

0.8 

2 

600 

0.6 

27 

0.4 

29 

0.4 

29 

0.2 

7 

0.5 

2 

800 

0.7 

24 

0.4 

27 

0.4 

27 

0.3 

7 

0.3 

2 

1200 

0.8 

11 

1.2 

23 

1.2 

23 

0.0 

4 

0.0 

1 

Position 

Depth 


SD 


Total  SD 


H  I         J 

Total  SD   Total  SD   Total  SD   Total 


00 

1.0 

4 

0.6 

2 

0.7 

3 

0.5 

3 

2.2 

7 

100 

0.8 

4 

0.5 

2 

1.1 

3 

0.0 

1 

1.3 

7 

200 

0.8 

2 

0.0 

1 

1.4 

3 

0.0 

1 

1.3 

6 

300 

0.9 

2 

0.0 

1 

1.2 

3 

0.0 

1 

1.3 

6 

400 

0.8 

2 

0.0 

1 

0.6 

3 

0.0 

1 

1.7 

6 

600 

0.5 

2 

0.0 

1 

0.3 

3 

0.0 

1 

2.2 

5 

800 

0.3 

2 

0.0 

1 

0.1 

3 

0.0 

1 

2.8 

4 

1200 

0.0 

1 

0.0 

1 

0.2 

3 

0.0 

1 

3.4 

4 
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Figure  2.   Graph  of  normalized  standard  deviation  for  data 

collection  points.   Values  are  obtained  from  Lewit's 
(1972)  tables.   Figure  A  is  for  areas  having  small 
standard  deviation,  figure  B  is  for  areas  with  large 
standard  deviation.   Standard  deviation  given  in  C. 
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imposed  by  this  investigator,  points  A,  B,  C  and  J  are  judged  as 
having  a  large  vertical  temperature  variability  and  while  at  points 
D,  E,  F,  G,  H  and  I,  the  temperature  variability  is  expected  to  be 
low.   As  expected,  points  near  strong  currents,  e.g.,  the  Gulf 
Stream,  show  large  variability  while  points  in  relatively  calm 
waters  e.g.,  the  Sargasso  Sea,  show  low  variability. 

Table  II  shows  the  comparison  of  the  FNWC  climatology 
(FNCL)  and  NAVOCEANO  climatology  (NAVOCL)  with  the  actual  mean 
temperature  profile  for  the  three  areas  A,  B,  and  J.   Differences 
between  each  climatology  and  the  mean  for  each  12-hour  forecasting 
period  are  given  along  with  the  normalized  figures.   Normalization 
is  accomplished  by  dividing  the  difference  at  any  level  by  the 
maximum  difference  of  either  column. 

Normalization  is  computed  to  establish  a  percentage  compari- 
son between  climatologies  of  FNWC  and  NAVOCEANO  with  respect  to  the 
mean  temperature  profile.   Figures  3  a,b,c  and  d  are  graphs  of  the 
normalized  differences  with  linear  interpolation  between  points  for 
positions  B,  A,  and  J.   Were  a  climatology  to  agree  perfectly  with 
the  mean,  the  curves  would  rest  exactly  on  the  zero  ordinate.   The 
relative  size  of  the  normalized  difference  at  any  level  is  a  measure 
of  the  error  associated  with  that  climatology  such  that  1.0  is  the 
maximum  relative  error.   Thus  an  area  between  the  curve  and  the  0.0 
normalized  difference  line  is  an  indication  of  the  error  present. 

AREA  FNCL  -  AREA  NAVOCL  =  DECREASE  IN  ERROR   (1) 

AREA  FNCL 

Equation  (1)  is  used  to  give  a  quantitative  measure  of  the  decrease 

in  error  associated  with  the  FNCL. 
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A 


Table  II.   Comparison  of  FNCL  and  NAVOCL  with  mean  temperature  pro- 
files.  Differences  are  normalised.   Number  of  observations 
upon  which  the  mean  is  calculated  is  given  in  parenthesis. 
Mean  profiles  are  given  for  indicated  Zulu  time  periods. 


Posit 

ion  A 

Depth 

FNCL  NAV0CT 

MEAN 

racL- 

NAVOCL 

(ft) 

(°c)  CO 

(4) 

Absolute 

Normal 

Absolute 

Normal 

0600 

0600 

0600 

0b00    ' 

0600 

00 

24.8  22.6 

24.7 

.2 

.1 

2.1 

1.0 

100 

24.3  21.8 

23.1 

.8 

.4 

1.3 

.6 

200 

20.7  20.6 

21.7 

1.0 

.5 

1.1 

.6 

300 

19.7  19.9 

20.7 

1.0 

.5 

.8 

.4 

400 

19.1  19.2 

19.8 

.7 

.3 

.6 

.3 

600 

18.6  18.2 

18.9 

.3 

.1 

.7 

.3 

800 

18.3  18.0 

18.5 

.2 

.1 

.5 

.1 

1200 

17.8  17.5 

17.6 

.1 

.0 

.1 

.1 

Position  J 


Depth 

FNCL 

NAVOCL 

MEAH 

FNCL-MEAN 

NAVOCL-MEAN 

(ft) 

(°C) 

(*C) 

(2) 

Absolute 

Normal 

Absolute 

Norma] 

1112 

1112 

1112 

1112 

1112 

' 0TT 

24.3 

22.5 

23.3 

1.0 

.5 

.8 

.4 

100 

24.3 

21.7 

22.3 

2.0 

1.0 

.6 

.3 

200 

20.5 

19.9 

19.8) 

.7 

.4 

.1 

.1 

300 

19.4 

18.8 

18.9 

.5 

.2 

.1 

.1 

400 

18.9 

18.8 

18.8 

.1 

.1 

.0 

.0 

600 

18.4 

18.0 

18.0 

.4 

.2 

.0 

.0 

800 

18.0 

17.7 

17.7 

.3 

.2 

.0 

.0 

1200 

17.5 

17.5 

17.2 

.3 

.2 

.3 

.2 

Position  B 

)epth 

FNCL 

NAVOCL 

MEAN 

FNCL-MEAN 

NAVOCL-MEAN 

(ft) 

(°C) 

(°C) 

(3) 

(3) 

Absolute 

Normal 

Absolute 

Normal 

0800 

0812 

0800  0812 

0800 

0812 

0800  0812 

0800  12 

00 

24.2 

22.5 

25.1 

25.1 

.9   .9 

.4 

.4 

2.6  2.6 

1.0  1 

.0 

100 

24.0 

22.5 

23.0 

23.6 

1.0  .4 

.4 

.1 

.5  1.1 

.2 

A 

200 

20.5 

20.0 

21.3 

21.6 

.8  1.1 

..3 

.4 

1.3  1.6 

.5 

( 

300 

19.5 

18.8 

20.7 

20.5 

1.2  1.0 

.5 

.4 

.9  1.7 

.4 

.1 

400 

18.9 

18.5 

20.2 

20.0 

1.3  1.1 

.5 

.4 

1.7   1.5 

.7 

.6 

600 

18.6 

18.0 

19.5 

19.2 

.9   .6 

.4 

.2 

1.5   1.2 

.6 

.5 

800 

18.6 

17.8 

18.9 

18.5 

.5   .1 

.2 

.0 

1.1    .7 

.4 

3 

^200 

17.9 

17.5 

18.3 

17.5 

.4   .4 

.1 

.2 

.8   .0 

•  3 

0 

Position  B 

Depth 

FNCL 

NAVOCL 

MEAN 

FNCL-MEAN 

NAVOCL-MEAN 

(ft) 

(°C) 

(°C) 

(3) 

(2) 

Absolute  Normal 

Absoulte   Normal 

0000 

0912 

0900  0912 

0900 

oyi2 

0900  0912 

0900  0912 

00 

24.2 

T2.5 

25.0 

25.0 

.8   .8 

.3 

•: 

2.5  2.5 

1.0  .8 

100 

24.0 

22.5 

22.9 

24.2 

1.1   .2 

.4 

.1 

.4  1.7 

.2   .6 

200 

20.5 

20.0 

20.2 

23.1 

.2   2.6 

.1. 

.9 

.7  3.1 

.3  1.0 

300 

19.5 

18.8 

20.5 

22.0 

1.0   .5 

.4 

.2 

1.7  1.2 

.7   .4 

400 

18.9 

18.5 

19.6 

20.3 

.7   .4 

.3 

.1 

1.1   1.8 

.4   .6 

600 

18.6 

18.0 

19.2 

19.1 

.6   .5 

.2 

.2 

1.2  1.1 

.5   .4 

800 

18.6 

17.8 

18.7 

18.8 

.3   .4 

.1 

.1 

.9  1.0 

.4   .3 

L200  |17.9 

17.5 

17.8 

17.8 

.1   .1 

.0 

.0 

.3   .3 

.1   .1 
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Figure  3A.   A  graph  of  the  normalized  difference  between  mean 

temperature  profile  and  NAVOCL  and  FNCL  for  position 
B  at  0800  local  time. 
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Figure  3B.   A  graph  of  the  normalized  difference  between  mean 

temperature  profile  and  NAVOCL  and  FNCL  for  position 
B  at  0900  local  time. 
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Figure  3C.   A  graph  of  the  normalized  difference  between  mean 

temperature  profile  and  FNCL  and  NAVOCL  for  position  A. 
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Figure  3D.   A  graph  of  the  normalized  difference  between  mean 

temperature  profile  and  FNCL  and  NAVOCL  for  position  J, 
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With  this  comparison  basis,  at  position  A  (Figure  3c)  FNCL 
yields  a  14%  smaller  error  than  NAVOCL.   Contrarily,  for  position  J 
(Figure  3d),  NAVOCL  yields  a  35%  error  reduction.   For  all  four 
forecast  periods  for  position  B  (Figure  3a, b)  FNCL  shows  a  substantial 
decrease  in  error  compared  to  NAVOCL.   Specifically,  it  yields  a 
29%,  52%,  51%  and  G0%  decrease  in  error  respectively. 

For  position  J  only  two  observations  are  available  to  obtain 
the  mean  temperature  profile.   For  the  other  positions  the  means  are 
based  upon  a  greater  number  of  observations.   Thus  for  position  J, 
that  NAVOCL  is  closer  to  the  mean  profile  is  probably  not  significant. 
For  all  remaining  forecasts,  FNCL  had  substantially  less  error  than 
NAVOCL.   It  is  concluded  from  this  trend  that  in  areas  of  expected 
large  standard  deviation,  FNCL  is  a  better  approximation  to  the  mean 
temperature  profile. 

Analysis  of  the  profiles  for  positions  D,  E,  G,  H  and  I 
(previously  found  to  be  in  low  standard  deviation  areas)  are  listed 
in  Table  III.   As  before,  the  differences  between  the  individual 
climatologies  and  the  mean  temperature  profile  are  normalized.   These 
normalized  values  are  plotted  in  Figures  4a-e  for  each  position. 

An  area  analysis  for  error  is  xjerf°rmed  as  for  high  standard 
deviation  positions.   In  the  low  standard  deviation  areas,  FNCL 
showed  less  error  than  NAVOCL  in  every  case  except  position  I. 
For  position  I,  NAVOCL  had  a  38%  smaller  error  than  FNCL,  the  reason 
for  which  is  unknown  to  this  investigator.   For  positions  D,  E,  G 
and  H  the  results  are  listed  below.   Percentages  given  are  the 
percent  of  decreased  error  experienced  by  FNCL  over  NAVOCL. 
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Table  III.   Comparison  of  FNWC  and  NAVOCL  with  mean  for  the  zulu 
time  period  indicated.   Differences  are  normalized. 
Number  of  observat ions  upon  wh ich  the  mean  is  calcu- 
lated is  given  in  parenthesis. 


Position  D 

Deptl 

FNCL 

NAV0CL 

MEAN 

FNCL-KEAN 

NAVOCL-MEAN 

(ft) 

CC) 

(°C) 

(6)        (3) 

Abso 

Lute        Normal 

Absoulte 

Morma 

I 

0900   U9T2 

0900 

091 J 

0900 

09  IS 

0900  0912 

0900 

0912 

00 

25.5 

22.5 

24.5   24.4 

1.0 

1.1 

.5 

.6 

2.0      1.9 

1.0 

1.0 

100 

24.8 

21.8 

22.9   23.2 

1.9 

1.6 

1.0 

.8 

1.1      1.4 

.6 

.7 

200 

21.3 

20.3 

21.3    21.3 

.0 

.0 

.0 

.0 

1.0      1.0 

.5 

.5 

300 

20.2 

18.8 

20.2    20.2 

.0 

.0 

.0 

.0 

1.4      1.4 

.7 

.7 

400 

19.6 

18.6 

19.7    19.5 

.1 

.1 

.0 

.0 

1.1        .9 

.6 

.5 

600 

18.9 

18.0 

19.0    18.9 

.1 

.0 

.0 

.0 

1.0        .9 

.6 

.5 

800 

18.4 

17.8 

18.6   18.5 

.2 

.1 

.1 

.0 

.8        .7 

.4 

•4 

1200 

17.7 

17.5 

17.7    17.8 

.0 

.1 

.0 

.0 

.2         .3 

.1 

•2    1 

Position  E 


Depth 

KNCL 

NAVOCL 

MEAN 

FNCL- MEAN 

NAVOCL- 

^EAN 

(ft) 

CC) 

CO 

(4) 

(8) 

Absolute 

Normal 

Absolute 

Normal 

1000 

1012 

1000   1012 

1000   12 

1000  1012 

1 000 

1012 

00 

25.  y 

^3.9 

24.2 

24.1 

1.7      1.7 

.5    .5 

.3      .3 

.1 

.1 

100 

25.9 

22.5 

22.4 

22.8 

3.5      3.1 

1  .  01 . 0 

.1       .3 

.0 

.1 

200 

21.9 

21.4 

21.1 

21.7 

.8         .2 

.2    .0 

.3      .3 

.1 

.1 

300 

20.6 

20.4 

20.0 

20.8 

.6         .2 

.1    .0 

.4      .4 

.1 

.1 

400 

19.8 

20.0 

19.7 

20.0 

.1         .2 

.0    .0 

.3      .0 

.1 

.0 

600 

18.9 

18.9 

18.7 

18.7 

.2         .2 

.0    .0 

.2      .2 

.0 

.0 

800 

18.3 

18.0 

18.1 

18.2 

.2         .1 

.0    .0 

.1      .2 

.0 

.0 

1200 

17.6 

17.5 

17.4 

17.4 

.2         .2 

.0    .0 

.1      .1 

.0 

.0 

Posit 

on   G 

Depth 

FNCL 

NAVOCL 

MEAN 

FNCL 

-MEAN 

NAVOCL-MF.AN 

(ft) 

CC) 

CC) 

(18) 

(9) 

Abso 

lute 

Normal     A 

>solute        Norma' 

TTTT 

1100 

1112 

1100 

1112 

1100 

1112 

1100   1112 

1100 

00 

25.4 

21.4 

24.3 

23.9 

1.1 

1.5 

.4 

.4 

1.9   2.5 

.8 

.7 

100 

25.4 

20.2 

22.4 

21.9 

3.0 

3.5 

1.0 

1.0 

2.2   1.7 

.7 

.5 

200 

21.6 

19.7 

20.9 

20.6 

.7 

1.0 

.2 

.3 

1.2      .9 

.4 

.3 

300 

20.4 

18.8 

20.2 

19.8 

.2 

.6 

.1 

.2 

1.4    1.0 

.5 

.3 

400 

19.7 

18.5 

19.7 

19.2 

.0 

.5 

.0 

.1 

1.2      .7 

.4 

.2 

600 

18.6 

18.0 

18.8 

18.2 

.2 

.4 

.  i 

.1 

.8      .5 

.3 

.2 

800 

18.0 

17.9 

18.3 

17.9 

.3 

.1 

.1 

.0 

.4      .0 

.1 

.0 

1200 

17.3 

17.3 

17.5 

17,1 

.2 

.2 

.1 

.0 

.2      .2 

.0 

.0 

Position  H 

Depth 

FNCL 

NAVOCI 

1*EAN 

FNCL-MEAN 

NAVOCL-MEAN 

(ft) 

("C) 

CO 

(2)      (9) 

Abso 

lute 

Normal 

Absolute      Normal 

1200  1212 

1200 

1212 

1200  121 

1200  1212 

1200  1212 

uu 

24.0 

20.0 

23. B  25.4 

.2 

1.4 

.0 

.2 

3.8      5.4 

1.0     1.0 

100 

24.0 

19.0 

22.2   22.8 

1.8 

1.2 

.5 

.2 

3.2      3.8 
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Figure  4A.   Normalized  difference  between  FNCL  and  NAVOCL  and 
the  mean  temperature  profile  for  position  D. 
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Figure  4B.   Normalized  differences  between  FNCL  and  NAVOCL 
and  the  mean  temperature  profile  for  position  E, 
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Figure  4C.   Normalized  difference  between  the  mean  temperature 
profile  and  the  NAVOCL  and  FNCL  for  position  G. 
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Figure  4D.   Normalized  difference  between  mean  temperature  profile 
and  the  FNCL  and  NAVOCL  for  position  H. 
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Figure  4E.   Normalized  difference  between  mean  temperature 
profile -and  the  NAVOCL  and  FNCL  for  position  I, 
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As  with  the  high  standard  deviation  areas,  FNCL  is  considered  a 
better  approximation  to  the  mean  temperature  profile  than  NAVOCL 
based  upon  the  strong  trend  observed. 
2.   Acoustic  Predictions 

Comparisons  of  temperature  profiles  and  related  conclusions 
about  their  validity  is  insufficient  for  a  complete  analysis.   Since 
temperature  profiles  are  not  the  final  x^roduct,  analyses  of  the 
acoustic  prediction  based  on  each  profile  are  warranted. 

Several  methods  of  comparison  of  the  acoustic  prediction 
based  on  climatology  and  the  mean  BT  data  are  available.   The  para- 
meter used  by  the  fleet  for  any  such  analysis  is  range;  therefore, 
this  analysis  will  center  on  predicted  range.   The  standard  set  of 
input  data  for  use  with  the  propagation  loss  program  are: 

wave  height  =  one  foot, 
frequency    =  60  Hz, 
source  depth  =  60  feet, 
receiver  depth  =  60  feet. 

With  this  approach  all  remaining  parameters  (bottom  depth  and  bottom 
roughness)  are  held  constant  except  the  temperature  profile.   This 
allows  for  a  comparison  of  the  affects  that  various  separate  tem- 
perature profiles  have  on  acoustic  propagation. 

The  publication,  ASW  Oceanographic  Environmental  Services 
(NAVAIR  50-1G-24,  1973) explains  the  method  and  logic  behind  the  figure 
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of  merit  (FOM)  determination  from  a  propagation  loss  curve.   This 
procedure  is  followed  here  with  a  FOM  value  of  90  db/ybar  assumed 
for  all  cases.   Using  the  ordinate  value  of  90  db  as  a  guide,  the 
range  to  the  first  intersection  with  the  propagation  loss  curve  is 
called  the  direct  path  range.   Should  convergence  zone  influence 
be  strong  enough,  the  90  db  line  curve  will  intersect  the  propaga- 
tion loss  curve  at  least  twice  more.   The  range  between  the  second 
and  third  intersection  is  the  width  of  the  convergence  zone  or 
annulus,  and  the  range  from  the  source  to  the  leading  edge  of  the 
annulus  is  the  range  to  the  first  convergence  zone.   Figure  5  shows 
an  example  of  the  method  for  range  determination.   These  three 
values  are  of  tactical  importance  to  the  fleet;  thus,  they  are 
pertinent  analysis  parameters. 

Table  IV  is  constructed  for  positions  B  and  J,  two  positions 
in  areas  of  high  standard  deviation,  and  H  and  E,  two  positions  in 
areas  of  low  standard  deviation.   Ranges  for  the  above  three  para- 
meters are  given.   Clark  (FNWC,  personal  communication)  indicates 
that  due  to  a  computer  roundoff  of  range  information,  accuracy  only 
to  the  nearest  nautical  mile  can  be  expected.   No  greater  accuracy 
is  claimed  nor  possible  by  normal  computer  output  means.   Such  condi- 
tions force  the  use  of  a  band  of  ranges  vice  a  single  range.   Com- 
patibility of  either  climatology  acoustic  forecast  with  the  acoustic 
forecast  based  on  the  mean  temperature  profile  is  achieved  should 
any  portion  of  the  band  of  ranges  computed  from  a  climatology  coin- 
cide with  any  portion  of  the  band  of  ranges  computed  from  the  mean 
temperature  profile. 
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Acoustic  range  analysis  for  high  standard  deviation 
positions  (B  and  J)  and  low  standard  deviation  pro- 
sitions  (H  and  E).   (I)  indicates  incompatibility  and 
(C)  indicates  compatibility  with  ranges  as  computed 
from  the  mean  BT  profile. 


Position  B 

FOM  =  90  db 

Acoustic 
Prediction 
Based  Upon 

Forecast 

Time 

Period 

Direct 
Path  (nm) 

First 
CZ(nm) 

Width 
of 
Annulus (nm) 

MBAH 

FNCL 

NAVOCL 

MEAN 
FNCL 
NAVOCL 

0800 

osoo 

0800 

0812 
0812 
0812 

1    -   4 
21  -  23  (I) 
21  -  23  (1) 

21  -  23 

21  -  23  (C) 

21  -  23  (C) 

35  -  37 

35  -  37(C) 
34  -  36(C) 

34  -  36 

35  -  37(C) 
34  -  36(C) 

1  -  3 
1  -  3(C) 

1  -  3(C) 

1  -  3 
1  -  3(C) 
1  -  3(C) 

Position  J 


FOM  =  90  db 


Acoustic 
Prediction 
Based  Upon 

Forecast 

Time 

Period 

Direct 
Path  (nm) 

First 
CZ(nm) 

Width 
of 

Annulus (np1 

MEAN 

rNCL 
NAVOCL 

MEAN 
FNCL 

NAVOCL 

1100 
1100 
1100 

1112 
1112 
1112 

16  -  18 

5   -   7  (1) 

2   -   4  (I) 

2   -   4  ■ 
2   -   4  (C) 
2   -   4  (C) 

No  CZ 

34  -  36(1) 

36  -  38(1) 

34  -  36 
34  -  36(C) 
36  -  38(C) 

nA 

1  -  3(1) 

1  -  3(1) 

1  -  3 
1  -  3(C) 

1  -  3(C) 

Pos  it  ion  H 

FOM  =  90  db 

Acoustic 
Prediction 
Based  Upon 

Forecast 

Time 

Teriod 

Direct 
Path  (nm) 

First 
CZ(nm) 

Width 

of 
\nnulus (nm' 

MEAN 

FNCL 
NAVOCL 

MEAN 
FNCL 
NAVOCL 

1200 
1200 
1200 

1212 
1212 
1212 

0 

5  -  7  (1) 

0  (C) 

0 

0  (C) 

0  (C) 

34  -  36 
34  -  36(C) 
37  -  39(1) 

34  -  36 
34  -  36(C) 
37  -  39(1) 

1  -  3 
1  -  3(C) 
1  -  3(C) 

1  -  3 
1  -  3(C) 
1  -  3(C) 

Posit  ion  E 

FOM  =  90  db 

Acoustic 
Fred  it  ion 
Based  Upon 

.  Forecast 
Time 
Period 

Direct 
Path  (nm) 

First 
CZ(nm) 

Width 
of 
Annulus (nm) 

MEAN 
FNCL 
NAVOCL 

MEAN 
FNCL 
NAVOCL 

1000 
1000 
1000 

1012 
1012 
1012 

2  -  4 

4  -  6  (C) 

2  -  4  (C) 

2-4 

4  -  6  (C) 

2  -  4  (C) 

34  -  36 
34  -  36(C) 
34  -  36(C) 

34  -  36 
34  -  36(C) 
34  -  36(C) 

no          on 

36 


a 
u 

3 

60 
•H 


PPQ   iJ  O  ifl  w 


37 


Despite  the  relative  superiority  of  FNWC  temperature  pro- 
files for  positions  B,  H  and  E,  it  is  necessary  to  analyze  the 
three  range  predictions  separately  since  different  influences 
affected  the  development  of  each.   In  this  way,  errors  inherent 
in  each  system  can  be  examined. 

For  position  H,  a  low  standard  deviation  area,  the  direct 
path  range  from  NAVOCL  is  compatible  with  the  mean,  but  not  so  for 
the  FNCL  direct  path  range.   At  position  E  compatibility  exists 
at  all  ranges  for  both  climatologies.   The  high  standard  deviation 
positions,  B  and  J,  show  compatible  direct  path  ranges  in  the 
final  forecasts,  but  are  incompatible  in  the  initial  forecast. 
3 .   Discussion  of  Results  \y 

One  of  the  main  contributers  to  the  magnitude  of  the 
direct  path  range  is  the  thickness  of  the  MLD.   For  example,  Harvey 
(1972)  shows  that  at  1700  Hz  the  energy  loss  for  a  MLD  of  100  feet 
is  2.3  db/nm  whereas  for  a  MLD  of  200  feet  the  loss  is  1.3  db/nm. 
This  MLD  thickness  effect  is  significant  since  over  a  typical  direct 
path  range  of  5  run,  5  db  difference  occurs.   Thus  an  error  in  MLD 
thickness  of  a  few  tens  of  feet  can  appreciably  alter  a  forecast. 
Also,  the  MLD  is  a  function  of  space  and  time.   Figure  6  shows  the 
MLD  from  analysis  of  BT  data  taken  in  SEACONEX.   The  crosshatched 
area  depicts  the  MLD.   Both  time  and  range  effects  are  present. 
The  cross-section  was  composed  of  BT  data  collected  at  different 
times  and  places.   The  horizontal  distance  was  32  run  and  a  time 
of  two  hours  to  travel  the  distance.   The  variability  of  the  thick- 
ness of  the  MLD  is  clearly  shown.   In  places,  a  change  of  50  feet 
in  the  MLD  can  occur  over  a  distance  of  only  a  single  nautical  mile. 
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Because  of  its  larger  data  base,  the  FNCL  was  expected  to 
show  a  better  correlation  to  the  actual  mean  MLD  than  did  NAVOCL. 
In  general,  the  temperature  profile  analysis  supports  this.   Yet, 
since  direct  path  range  is  a  strong  function  of  MLD  (Urick,  1967), 
the  acoustic  forecast  indicates  no  such  superiority.   Neither  cli- 
matology showed  any  superiority  in  acoustic  compatibility.   For 
positions  B,  J,  and  E  almost  identical  records  of  compatibility 
are  achieved.   Only  for  position  H  did  a  discrepancy  exist. 

Any  climatology  is,  in  essence,  a  mean  with  an  associated 
standard  deviation.   The  actual  real-time  MLD  could  lie  anywhere 
within  certain  limits  about  that  mean  MLD.   The  size  of  the  limits 
is  dependent  upon  the  magnitude  of  the  standard  deviation.   For 
example,  should  the  distribution  of  MLD  be  Gaussian,  there  would 
be  a  68%  probability  that  the  actual  MLD  would  lie  within  one 
standard  deviation  of  the  mean  MLD.   The  actual  MLD  then  will  lie 
within  an  envelope  about  the  mean  MLD. 

For  a  short  time  period,  e.g.,  one  or  two  days,  the  actual 
MLD  can  lie  anywhere  within  that  envelope.   As  long  as  the  FNCL  and 
the  NAVOCL  acoustic  predictions  lie  within  that  envelope,  they  have 
•equal  chance  of  correlation  to  the  actual  MLD.   Over  a  long  time 
period  the  actual  MLD  will  vary  about  the  mean  MLD.   It  is  the 
opinion  of  this  investigator  that  even  though  a  particular  climato- 
logy over  a  long  time  period  may  provide  a  better  representation  of 
the  actual  MLD  than  another  climatology,  over  any  short  time  period 
no  superiority  is  observed.   Direct  ranges  predicted  by  any  climato- 
logy must  be  considered  as  a  "first  guess"  until  additional  data 
is  available. 
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The  range  to  the  first  convergence  zone  and  width  of  the 
annulus  is  a  function  of  the  shape  of  the  SVP  and  the  magnitude 
of  the  gradient  of  the  thermocline.   Again  it  is  expected  that  the 
larger  data  base  and  more  sophisticated  thermal  structure  of  FNCL 
would  allow  the  ranges  associated  with  it  to  be  closer  than  the 
ranges  associated  with  NAVOCL  when  compared  with  ranges  obtained 
by  the  actual  mean  temperature  profile.   As  with  the  direct  path 
case,  no  such  superiority  was  seen  in  the  acoustic  range  predic- 
tions.  For  these  ranges,  the  FNCL  showed  excellent  compatibility 
with  the  mean  at  all  positions  except  the  initial  forecast  of 
position  J.   Even  for  this  case,  the  ranges  forecasted  were  no 
farther  from  the  mean  forecast  than  those  of  NAVOCL.   For  the 
NAVOCL,  incompatibility  is  experienced  at  positions  H  and  J. 

Again  this  investigator  considers  that  conditions  on  a 
daily  basis  are  too  variable  for  either  climatology  SVP  to  yield 
consistently  better  results  than  the  other.   Local  effects  could 
influence  the  mean  temperature  profile  significantly  from  the 
climatological  mean.   On  the  average  the  results  of  any  correla- 
tion between  any  two  ranges  is  considered  only  a  temporary  effect 
and  not  to  be  expected  for  subsequent  reports.   This  is  seen  for 
position  J,  which  has  no  correlation  for  FNCL  in  the  first  fore- 
cast, but  exact  matching  in  the  second.   Figure  6  shows  a  tongue 
of  cold  water  that  extends  from  depth  to  near  surface.   According 
to  Nan-Titi  (1967) ,  this  tongue  may  be  only  a  local,  short-lived 
effect  (a  few  days)  that  would  be  reflected  in  a  local  mean  but  not 
in  the  climatology.   This  is  further  evidence  that  a  monthly  mean 
cannot  be  expected  to  correlate  exactly  with  daily  observations. 
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This  investigator  believes  that  for  the  purpose  of  daily- 
forecasting  no  advantage  is  gained  from  the  sophisticated  approach 
to  climatology  used  by  FNWC  for  the  case  of  a  transiting  ship. 
Should  a  lengthy  time  (30  days  or  more)  be  spent  in  a  single  area, 
the  FNCL  forecast  may  prove  to  give  better  results  than  NAVOCL. 
When  only  a  few  days  are  spent  in  an  area,  the  daily  fluctuations 
about  the  average  will  make  correlations  with  any  climatology  on]y 
a  chance  effect. 

C.   SYNOPTIC  COMPARISON 

1.   FNWC  Synoptic  Temperature  Profile 

Every  BT  taken  during  "SEACONEX"  was  sent  to  FNWC,  where 
it  was  then  included  in  the  mean  synoptic  temperature  profile  as 
explained  in  Chapter  I.   Propagation  loss  curves  were  run  for  the 
positions  of  interest  (same  positions  as  used  in  previous  climato- 
logy analysis)  for  each  12-hour  analysis  period. 

The  FNWC  thermal  analysis  forecast  at  0000Z  is  intended  to 
be  valid  during  the  period  0000Z-1200Z.   This  profile  was  compared 
to  the  mean  of  all  BT's  taken  during  the  interval  0000Z-1200Z. 
This  mean,  constructed  of  on-scene  BT's,  was  considered  to  be  the 
actual  mean  for  the  period.   The  same  procedure  is  used  for  the 
1200Z  forecast.   No  control  is  placed  on  the  generation  of  the  mean 
profile  in  that  all  BT's  are  included  and.no  weighting  or  biasing 
is  used. 

Table  V  lists  the  temperature  profile  as  analyzed  by  FNWC 
for  the  time  periods  indicated.   The  notation  in  the  table  0600Z 
should  be  interpreted  as  0000Z  on  the  sixth  of  June.   A  major 
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Table  V.   Temperature  (C)  at  various  levels  for  FNWC  mean  temper- 
ature profile  and  the  actuel  mean  temperature  profile. 
The  difference  between  the  two  are  also  given.   Zulu 
time  periods  are  given. 


Posit  ion  A 


Depth 

FNWC 

Actual 

Difference 

(ft) 

mean 

mean 

FN-Actual 

06002 

0600Z 

0600Z 

00 

24.8 

24.7 

0.1 

100 

21.0 

23.1 

-2.1 

200 

19.8 

21.7 

-1.9 

300 

18.8 

20.7 

-1.9 

400 

18.6 

19.8 

-1.2 

600 

18.1 

18.9 

-0.8 

800 

17.6 

18.5 

-0.9 

1200 

17.3 

18.6 

-0.3 

Pi'Fiti 

>n  B 

Depth 

FNWC 

Actual 

Difference 

(ft) 

mean 

mean 

FN-Actual 

0800 

0812    0900   0412 

0800   0812   0900 

0912 

0800   0812 

0900   0912 

25.0   25.0   24.8 

25.1    25.1    25.0 

25.1 

-0.6   -0.1 

-0.0  -0.2 

100 

24.3 

24.0  23.0  24.0 

23.0  23.6   22.9 

24.21 

-1.3   -0.4 

-0.1    -0.2 

200 

21.8 

20.7  20.9  20.5 

21.6   21.3   20.7 

23.1 

0.2    -0.6 

0.2   -2.6 

300 

20.2 

19.7   19.9    19.5 

20.7   20.5    20.5 

22.0 

-0.5    -0.8 

-0.6    -2.5 

400 

19.9 

19.1   19.6  19.2 

20.2  20.0  19.6 

20.3 

-0.3   -0.9 

-0.0   -1.1 

600 

19.3 

18.5   19.0   18.5 

19.5    19.2    19.2 

19.1 

-0.2    -0.7 

-0.2    -0.6 

800 

19.1 

18.0  18.5   18.0 

18.9   18.5   18.7 

18.8 

0.2    -0.5 

-0.2    -0.8 

1200 

18.8 

17.0  17.7    17.6 

18.3    17.5   17.8 

17. fi 

0.5    -0.5 

-0.1    -0.2 

Position  F. 


Position  D 


Depth 

FNWC 

Actual 

Difference 

(ft) 

mean 

mean 

FN-Actual 

1000Z    1012Z 

1000Z 

10127. 

1000Z 

1012Z 

uu 

25.0     24.7 

24.2 

24.2 

0.8 

'  0.5  ' 

100 

22.8      23.2 

22.4 

22.8 

0.4 

0.4 

200 

21.4      21.7 

21.1 

21.7 

0.3 

0.0 

300 

21.0     21.2 

20.0 

20.8 

1.0 

0.4 

400 

20.6     20.5 

19.7 

20.0 

0.9 

0.5 

600 

20.3     19.1 

18.7 

18.7 

1.6 

0.4 

800 

19.2     18.1 

18.1 

18.2 

1.1 

-0.1 

1200 

18.3     17.8 

17.4 

17.4 

0.9 

0.4. 

Depth 

FNWC 

Actual 

difference 

(ft) 

mean 

mean 

FN- 

Actual 

0900   0912 

0900   0912 

0900 

0912 

00 

25.5   24.8 

24.5    24.4 

1.0 

0.4 

100 

23.7   23.8 

22.9   23.2 

0.8 

0.6 

200 

22.3   21.7 

21.3   21.3 

1.0 

0.4 

300 

21.4   20.9 

20.2    20.2 

1.2 

0.7 

400 

20.6   20.4 

19.7   19.5 

0.9 

0.9 

600 

19.7  19.5 

19.0   18.9 

0.7 

1.6 

800 

19.1   18.7 

18.6   18.5 

0.5 

0.2 

1200 

17.9   17.8 

17.7   17.8 

0.2 

0.0 

Position  C 


Depth 

FNWC 

Actual 

Dlf  ferenct! 

(ft) 

mean 

mean 

FN-Aclual 

1100Z    11127 

11007. 

1112Z 

11G0Z   11122 

00 

24.6  23.6 

24.3 

23.9 

'    0.3      -0.3 

100 

23.3  22.9 

22.4 

21.9 

0.9        1.0 

200 

21.2    21.0 

20.9 

20.6 

0.3        0.4 

300 

20.3   20.1 

20.2 

19.8 

0.1        0.3 

400 

19.5  19.4 

19.7 

19.2 

-0.2        0.2 

600 

18.0  18.5 

18.8 

18.2 

-0.8        0.3 

800 

17.3  17.9 

18.3 

17.9 

-1.0        0.0 

UW 

16.6  16.7 

17.5  _ 

17.1 

-0.9        -0.4 

Positl 

Dn  H 

Depth 

FNWC 

Actual 

Difference 

(ft) 

mean 

mean 

FN-Actual 

1200Z    12127 

1200Z 

1212Z 

12O0Z    1212Z 

UU 

23.1      23.3 

23.8 

25.4 

-0.7      -2.1 

100 

22.8      21.9 

22.2 

22.8 

0.6      -0.9 

200 

22.0     20.5 

20.6 

21.1 

1.4      -0.6 

300 

20.8      19.7 

19.9 

20.6    • 

0.9      -0.9 

400 

19.2     19.1 

19.4 

20.2 

-0.2      -0.9 

600 

18.0     18.5 

18.3 

18.7 

-0.3      -0.2 

800 

17.4      18.1 

17.6 

17.8 

-0.2        0.3 
-0.3        0.1     i 

1200 

16.3     17.3 

16.6 

17.2 

Positi 

TO   I 

Depth 

FNWC 

Actual 

Difference 

(ft) 

mean 

mean 

FN-Actual 

1300Z 

1300Z 

1300Z 

00    " 

24.9 

26.0 

'  -1.1 

100 

23.8 

24.7 

-0.9 

200 

21.7 

22.5 

-0.8 

300 

21.1 

20.8 

0.3 

400 

20.4 

19.8 

0.6 

600 

18.8 

18.3 

0.8 

800 

17.7 

17.9 

-0.2 

1200 

16.7 

17.1    .. 

-0.4 

Position  J 


Dc  pth 

FNWC 

Actual 

Difference 

(ft) 

mean 

r:ean 

FN-Actual 

1112Z 

1112Z 

1112Z 

00 

23.4 

23.3 

-0.1 

100 

22.3 

22.3 

0.0 

200 

20.5 

19.8 

0.7 

300 

19.8 

18.9 

0.9 

400 

19.3 

18.8 

0.5 

600 

18.5 

18.5 

0.5 

800 

17.8 

17.7 

0.1 

1200 

16.3 

1   17.2 

-0.9 
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contribution  to  this  profile,  according  to  Chapter  II,  is  the  BT's 
taken  in  "SEACONEX."   However,  BT's  taken  during  the  0000Z-1200Z 
period  are  not  utilized  by  FNWC  until  the  1200Z  analysis.   Thus 
the  influence  of  BT's  taken  during  one  forecast  period  is  not  felt 
until  the  next  forecast  period.   Conditions  that  influenced  those 
BT's  may  not  be  present  during  the  next  forecast  period.   Hesse 
(1973)  suggested  that  ideally  as  more  BT  data  are  sent  to  FNWC  to 
be  included  in  the  mean  temperature  profile,  the  analyzed  mean 
should  better  approximate  the  actual  mean  temperature  profile.   On 
the  premise  that  temperature  changes  in  the  ocean  are  slow  (Nan-Titi, 
1967)  the  BT  data  sent  to  FNWC  could  be  used  despite  the  time  lapse 
between  taking  a  BT  and  its  utilization. 

One  way  to  test  the  validity  of  the  assumption  that  increas- 
ing the  BT  data  base  for  FNWC  analyses  will  force  the  temperature 
profile  to  converge  to  the  actual  temperature  profile,  is  to  compare 
two  actual  forecasts  for  the  same  position,,   Table  V  shows  the  differ- 
ence between  the  FNWC  mean  temperature  profile  and  the  actual  mean 
profile.   According  to  the  above  assumption,  this  difference  should 
decrease  with  each  succeeding  forecast.   Figure  7  is  a  plot  of  the 
difference  between  the  FNWC  mean  temperature  profile  and  the  actual 
mean  temperature  profile. 

Only  for  position  E  can  a  definite  pattern  be  seen  for  con- 
vergence at  all  eight  levels,  i.e.,  at  all  levels  the  difference 
values  for  the  later  analysis  were  smaller  than  those  for  the  initial 
analysis.   For  positions  B,  D,  G  and  H  46%  of  all  observations  showed 
a  greater  difference  in  the  later  analysis  as  compared  to  the  preced- 
ing one.   Conversely,  45%  have  a  lesser  difference  and  9%  experience 
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Figure  7A.   Differences  between  FNWC  mean  temperature  profile  and 
actual  mean  temperature  profile  are  plotted  for  the 
time  periods  indicated  at  position  B. 
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Figure  7B.  Difference  between  FNWC  mean  temperature  profile  and 
actual  mean  temperature  profile  are  plotted  for  time 
periods  indicated  at  position  D. 
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Figure  7C.   Differences  between  FNWC  mean  temperature  profile  and 
actual  mean  temperature  profile  are  plotted  for  the 
time  periods  indicated  at  position  G. 
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Figure  7D.   Differences  between  FNWC  mean  temperature  profile  and 
actual  mean  temperature  profile  are  plotted  for  the 
time  periods  indicated  at  position  H. 
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Figure  7E.   Differences  between  TNWC  mean  temperatui-e  profile  and 
actual  mean  temperature  profile  are  plotted  for  the 
time  periods  indicated  at  position  E. 
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no  change.   No  pattern  of  general  divergence  or  convergence  for 
any  level  can  be  established. 

Limited  data  do  not  allow  any  definitive  conclusions  as  to 
convergence  or  divergence.   The  data  that  are  available  show  that 
even  in  areas  of  small  temperature  variability,  the  assumption  that 
convergence  to  the  mean  temperature  profile  can  be  expected  with 
increasing  data  is  unfounded.   Apparently  the  problem  of  time  delay 
from  observation  to  data  input  by  FNWC  is  critical.   The  temperature 
profile  analysis  is  updated  by  FNWC  every  12  hours.   Therefore, 
any  BT's  taken  could  be  delayed  up  to  12  hours  or  more  before 
entry  into  the  analysis.   Small  scale,  local  anomalies,  which  may 
have  influenced  the  BT,  can  change  rapidly  making  the  BT  informa- 
tion erroneous.   This  is  compounded  by  the  problem  of  a  transiting 
ship  since  a  spatial  change  will  also  come  into  effect. 
2.   FNWC  Acoustic  Forecast 

Propagation  loss  profiles  for  positions  H,  E  and  B  were  used 
to  establish  a  range  analysis  based  on  the  FNWC  synoptic  mean  tem- 
perature profile  and  the  actual  mean  temperature  profile.   The  pro- 
cedure was  similar  to  that  used  in  Chapter  II.   According  to  Figure  7, 
at  position  B  the  temperature  profile  converged  for  the  upper  100  feet, 
but  was  divergent  below  100  feet.   Position  H  was  an  irregular  profile 
of  convergence  and  divergence.   An  entirely  convergent  pattern  was 
experienced  at  position  E.   Analysis  of  these  different  temperature 
patterns  will  permit  one  to  study  the  effects  that  different  profiles 
have  on  the  acoustic  ranges. 

Table  VI  shows  the  forecasted  acoustic  ranges  based  upon 
actual  mean  temperature  profiles  and  FNWC  mean  profiles.   All  FNWC 
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Table  VI.   Range  compatibility  between  acoustic  forecasts  based 
on  a  mean  synoptic  temperature  profile  and  the  FNWC 
mean  temperature  profile.   (C)  denotes  compatibility, 
(I)  denotes  incompatibility. 


Position  B 

FOM  = 

90  db 

Acoustic 
Prediction 
Based  Upon 

Forecast 

Time 

Period 

Direct 
Path(nm) 

First 
CZ(nm) 

Width 
of 
Annulus(nm) 

BT  Mean 
FNWC  Mean 

BT  Mean 
FNWC  Mean 

0800 
0800 

0812 
0812 

2-4 
21  -  23  (I) 

21  -  23 

21  -  23  (C) 

35  -  37 

35  -  37  (C) 

34  -  36 

34  -  36  (C) 

1-3 

1  -  3  (C) 

1  -  3 

1  -  3  (C) 

Position  H 

FOM  =  90  db 

Acoustic 

Forecast 

Direct 

First 

Width 

Prediction 

Time 

Path(nm) 

CZ(nm) 

of 

Based  Upon 

Period 

Annulus(nm) 

"BT  Mean 

1200 

0  -  1 

34  -  36 

1  -  3 

FNWC  Mean 

1200 

0  -  1  (C) 

34  -  36  (C) 

1  -  3  (C) 

BT  Mean 

1212 

0-1 

34  -  36 

1  -  3 

FNWC  Mean 

1212 

0  -  1  (C) 

35  -  38  (C) 

3  -  5  (C) 

Position  E 

FOM  -  90  db 

Acoustic 

Forecast 

Direct 

First 

Width 

Prediction 

Time 

Path(nm) 

CZ(nm) 

of 

Based  Upon 

Period 

Annulus  (nrr.) 

"BT  Mean 

1000 

2  -  4 

34  -  36 

I  -  3 

FNWC  Mean 

1000 

2  -  4  (C) 

34  -  36 

(C) 

1  -  3  (C) 

BT  Mean 

1012 

2  -  4 

34  -  36 

1-3 

FNWC  Mean 

1012 

2  -  4  (C) 

34  -  36 

(C) 

1  -  3  (C) 
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forecasts  were  based  on  at  least  one  BT  taken  in  "SEACONEX."   An 
analysis  similar  to  the  one  for  ranges  based  on  climatology  in  Chapter  II 
Section  B.2  was  used.   An  envelope  of  ranges  was  considered  for 
the  three  parameters  of  direct  path  range  to  the  first  convergence 
zone  and  width  of  the  convergence  zone  annulus.   By  definition 
compatibility  is  achieved  should  any  part  of  one  forecast  band 
overlap  with  the  other. 

Position  B  indicates  a  considerable  error  in  the  direct  path 
range  for  the  initial  forecast  period,  but  compatibility  is  achieved 
for  the  range  to  the  first  convergence  zone  and  width  of  the 
annulus.   By  the  time  of  the  second  forecast,  compatibility  is 
achieved  for  all  ranges.   Positions  H  and  E  show  complete  compati- 
bility at  all  times  and  ranges. 

These  developments  indicate  that  regardless  of  whether  tem- 
perature profiles  converge  or  diverge  to  a  mean,  once  FNWC  has 
received  several  (more  than  one)  BT's  from  an  area  the  ranges  pre- 
dicted are  generally  compatible  with  those  based  on  a  mean  synoptic 
profile.  Compatibility  cannot  be  assumed  for  a  single  BT  input  as 
seen  for  position  B  which  was  based  upon  only  a  single  BT  input. 
The  direct  path  range  was  incompatible  with  the  mean  direct  path 
range.   All  later  forecasts  were  based  upon  at  least  two  forecasts. 
Since  complete  compatibility  was  observed  for  these  later  forecasts, 
it  is  expected  that  ranges  predicted  by  FNWC  are  reliable  within 
roundoff  error  limits,  provided  at  least  two  BT's  have  been  input 
into  the  analysis. 
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III.   I CAPS  I/O  ANALYSIS 

A.   TEMPERATURE  INPUT 

ICAPS  operates  on  a  single  BT  input  upon  which  the  acoustic 
propagation  loss  profile  is  generated.   In  order  to  establish  how 
much  reliability  can  be  assigned  to  a  single  BT,  Table  VII  is 
presented.   Eor  the  construction  of  this  table,  BT's  taken  at  each 

position  are  averaged  and  a  variance  is  calculated  for  every  level. 

2 

Four  positions,  B,  D  ,  E  and  G,  have  a  variance  (a  )  greater  than 

1.0°C  for  at  least  one  level.  Positions  H  and  I  have  at  least  one 
level  which  has  a  variance  greater  than  0.5°C.  Thus  only  three  of 
nine  positions  show  a  variance  at  all  levels  less  than  0.5°C. 

The  above  analysis  indicates  that  in  all  geographic  positions 
examined,  a  high  degree  of  thermal  structure  variance  was  present. 
Thus  at  a  given  location  the  temperature  profile  can  vary  signifi- 
cantly over  a  relatively  short  period  of  time.  A  change  of  0.5  to 
1.0  degrees  centigrade  at  various  levels  can  alter  the  MLD  and  the 
gradient  of  the  thermocline  presenting  different  characteristics  for 
separate  profiles. 

Figure  6  shows  that  over  a  very  short  distance  of  the  ocean  the 
temperature  profile  can  vary  considerably.   For  the  profile  taken  at 
position  1  a  40  foot  MLD  and  a  moderately  steep  thermocline  are 
observed.   At  position  2  the  trace  shows  no  MLD  and  a  slightly  steeper 
thermocline.   The  separation  between  the  stations'  is  only  two 
nautical  miles  and  8  minutes  in  space  and  time. 
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Table  VII.   Variances  (aL)    calculated  for  all  positions  based 
on  actual  BT  data  taken  in  "SEACONEX". 


Position 

A 

Depth 
(ft) 

2 
a 

00 

.38 

100 

.17 

200 

.45 

300 

.14 

400 

.18 

600 

.05 

800 

.04 

1200 

.04 

Position 

B 

Depth 

2 

(ft) 

00 

.14 

100 

2.07 

200 

.44 

300 

.28 

400 

.38 

600 

.17 

800 

.17 

1200 

.17 

Position 

C 

Depth 
(ft) 

2 

a 

00 

.28 

100 

.17 

200 

.05 

300 

.04 

400 

.15 

600 

.00 

800 

.01 

1200 

.02 

Position 

D 

fcepth 

2 

(ft) 

a 

00 

.22 

100 

.60 

200 

1.08 

300 

.39 

400 

,26 

600 

.30 

800 

.21 

1200 

.20 

Position 

E 

Depth 

2 

(ft) 

a 

00 

.10 

100 

.96 

200 

1.08 

300 

.70 

400 

.68 

600 

.70 

800 

.69 

1200 

.06 

Position  G 

Depth 

2 

(ft) 

a 

00   . 

.39 

100 

1.95 

200 

1.15 

300 

.80 

400 

.52 

600 

.34 

800 

.12 

1200 

.26 

Position 

H 

Depth 

2 

a 

(ft) 

00 

.33 

100 

.40 

200 

.30 

300 

.53 

400 

.51 

600 

.07 

800 

.14 

1200 

.28 

Position 

I 

Depth 

2 
a 

(ft) 

00 

.11 

100 

.92 

200 

.17 

300 

.10 

400 

.22 

600 

.07 

BQ0 

.02 

1200 

.05 

Position 

J 

Depth 

2 

a 

(ft) 

00 

.02 

100 

.16 

200 

.02 

300 

.09 

400 

.06 

600 

.00 

800 

.02 

L200 

.06 
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Effects  of  internal  waves,  turbulence,  and  local  anomalies  can 
be  expected  to  affect  any  BT  profile.  The  mean  temperature  profile 
of  an  area  is  considered  as  the  best  representation  of  the  area,  but 
the  various  influences  present  can  limit  any  single  profile  in  its 
correlation  to  that  mean.   Figures  8a, b  and  c  are  constructed  to 
show  the  relationship  between  a  mean  profile  and  the  individual  BT's 
that  were  used  to  establish  it.   Positions  A  and  H  show  significant 
differences  between  BT's  with  respect  to  gradient  and  gradient 
changes.   For  the  most  part  there  is  little ' similarity  to  the  mean 
for  any  single  profile.   For  position  B,  curves  (4) ,  (2)  and  (3) 
do  exhibit  fairly  uniform  and  reasonable  similarity  to  one  another 
and  to  the  mean.   Below  400  feet  all  profiles  have  a  gradient  very 
close  to  0.3  °C/100  feet.   Even  for  this  case,  curve  (1)  has  a 
MLD  of  100  feet  while  the  mean  profile  does  not  exhibit  any  MLD. 

If  the  mean  temperature  profile  is  assumed  to  be  a  reliable  esti- 
mate for  the  true  mean  profile  in  an  area,  this  investigator  con- 
cludes that  a  single  BT  is  generally  an  unreliable  estimate  of  the 
mean  temperature  profile  for  an  area. 

B.   I  CAPS  ACOUSTIC  OUTPUT 

As  with  analysis  of  climatological  and  synoptic  information, 
temperature  analysis  alone  is  inconclusive.   The  acoustic  informa- 
tion generated  by  the  profiles,  being  the  end  product,  is  the  final 
rule  against  which  correlation  is  measured. 

Table  VIII  is  constructed  for  positions  A,  B,  and  H  of  which 
the  first  two  are  positions  of  high  standard  deviation  and  the 
latter  of  low  standard  deviation.   Using  the  same  analysis  for  range 
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Figure  8A.   Temperature  profiles  for  the  mean  and  single  BT's 
are  given  for  position  H. 
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Figure  8B.   Temperature  profiles  for  the  mean  and  single  BT's 
are  given  for  position  B. 
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Figure  8C.   Temperature  profiles  for  the  mean  and  single  BT' s 
are  given  for  position  A. 
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Table  VIII.   Single  BT  input  to  ICAPS  acoustic  ranges  versus 

the  acoustic  ranges  for  the  actual  mean  temperature 
profile.   (C)  denotes  compatibility,  (I)  denotes 
incompatibility. 


Position  B 

FOM  = 

9C 

db 

Acoustic 

Forecast 

Direct 

First 

Width 

Prediction 

Time 

Path(nm) 

CZ(nm) 

of 

Based  Upon 

Period 

Annulus (nm) 

ICAPS 

Single  BT 

0812 

1  -  3  (C) 

35  -  37  (C) 

4  -  6  (C) 

Mean 

0812 

1  -  3 

35  -  37 

4  -  6 

Position  A 

FOM  - 

90 

db 

Acoustic 
Prediction 
Based  Upon 

Forecast 

Time 

Period 

Direct 
Path  (nm) 

First 
CZ(nm) 

Width 
of 
Annulus (nm) 

ICAPS 
Single  BT 

'lean 

0600 
0600 

1  -  3  (C) 

2  -  4 

33  -  35  (C) 
34--  36 

1  -  3  (C) 
1-3 

Position  H 

FOM  = 

90 

db 

Acoustic 

Forecast 

Direct 

First 

Width 

Prediction 

T  ime 

Path(nm) 

CZ(nm) 

of 

Based  Upon 

Period 

Annulus  (nm 

ICAPS 

Single  BT 

1200 

1  -  3  (C) 

34  -  36  (C) 

1  -  3  (C) 

'lean 

1200 

0-2 

34  -  36 

1-3 
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as  used  previously,  it  is  seen  that  for  the  three  cases  generated, 
complete  compatibility  for  the  acoustic  ranges  is  experienced. 

This  investigator  concludes  that  at  least  for  the  cases  shown, 
a  single  BT  can  be  a  reliable  indicator  for  acoustic  ranges.   The 
variability  in  temperature  profiles  has  little  influence  over  the 
acoustic  transmission  loss. /  Clark  (FNWC,  personal  communication) 

indicated  that  low  frequency  signals  tend  to  be  little  effected  by 

] 

minor  changes  in  the  temperature  profile./'  In  other  words,  minor 

changes  in  the  profile  are  smoothed  out  when  a  low  frequency  signal 
is  transmitted.   Consequently,  the  influences  of  small  scale  per- 
turbations in  the  temperature  profile  will  not  be  realized  for  low 
frequencies.   As  long  as  profiles  have  the  same  general  shape, 
approximately  equal  detection  ranges  can  be  expected. 

With  such  a  small  statistical  base,  conclusions  covering  any 
large  scale  area  are  unfounded.   Yet  the  evidence  that  is  available 
does  indicate  that  for  the  most  part  in  any  single  area,  a  single 
BT  input  can  furnish  the  ICAPS  operator  with  x-eliable  range  data. 
The  ranges  calculated  can  at  least  be  used  as  a  good  first  approxi- 
mation to  those  of  the  mean  within  normal  error  limits. 
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IV.   CONCLUSIONS 


A.  CLIMATOLOGY 

The  data  presented  show  that  for  a  large  percentage  of  areas, 
the  FNWC  climatological  temperature  profile  appeared  to  be  closer 
to  the  mean  temperature  profile  for  the  studied  areas  than  the  cli- 
matology developed  and  used  by  NAVOCEANO  for  the  ICAPS  computer 
for  the  exercise  "SEACONEX."   It  should  be  remembered  that  this 
NAVOCEANO  climatology  was  only  meant  as  a  first  attempt  until  a 
more  sophisticated  model  could  be  developed.   However,  low  frequency 
acoustical  forecasts  based  upon  either  climatology  do  not  indicate 
a  clear  preference  for  one  system  over  the  other.   While  local  mean 
temperature  profiles  will  smooth,  but  still  incorporate,  local 
anomalies,  low  frequency  sound  will  not  respond  to  small  temperature 
fluctuations. 

This  investigator  concludes  that  for  a  short  time  period  the 
assumptions  that  FNWC  climatology  will  yield  more  reliable  acoustic 
ranges  than  NAVOCEANO  climatology  is  unfounded.   A  transiting  naval 
ship  can  expect  nearly  the  same  results  from  a  sophisticated  climato- 
logy as  well  as  from  a  simple  one.   Forecasts  based  on  climatology 
alone  should  be  used  only  as  an  initial  approximation  and  only  for 
information  purposes.   The  probability  that  the  predicted  ranges  will 
reflect  actual  conditions  are  uncertain. 

B.  SYNOPTIC 

v/  FNWC  synoptic  temperature  profiles  show  no  overall  evidence  of 
convergence  to  a  mean  temperature  profile.   Yet,  within  limits  of 
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roundoff  error,  the  ranges  forecast  by  FNWC  based  upon  at  least 
two  BT's  proved  to  be  similar  to  that  based  upon  the  ir.ean  BT  pro- 
file.  It  is  concluded  that  FNWC  synoptic  forecasts  do  indeed 
provide  a  reliable  tactical  forecast. 

The  single  BT  input  for  ICAPS  has  proven  to  be  a  reliable 
forecast  agent.   However,  the  user  must  attempt  to  calculate  a 
mean  temperature  profile  for  his  area  as  soon  as  possible.   Until 
a  mean  profile  is  established,  a  single  BT  input  may  be  used  with  . 
moderate  reliability. 

C.   OPERATIONAL  CONSIDERATIONS 

Based  on  information  from  "SEACONEX,"  the  large  time  delay  (up 
to  12  hours)  between  the  time  a  BT  was  taken  and  incorporation  in 
the  temperature  profile  at  FNWC  and  the  time  delay  (average  of  two 
hours)  between  the  time  the  FNWC  forecast  was  made  and  receipt  of 
the  forecast,  and  the  fact  that  in  any  combat  situation,  communica- 
tions with  FNWC  will  be  minimal,  the  dependence  upon  FNWC  to  provide 
real-time  forecasts  to  the  fleet  seems  unrealistic.   ICAPS  utilizes 
the  same  acoustic  models  and  inputs  without  the  communications 
problems  with  a  very  short  time  delay.   It  is  concluded  that  the 
ICAPS  system  is  the  appropriate  real-time  system,  but  FNWC  can  be 
used  as  an  accurate  backup  system. 

For  tactical  planning  purposes,  FNWC  climatology  should  be  used 
as  the  source  of  acoustical  range  information,,   Their  large  choice 
of  data  presentation  techniques  and  ability  to  widely  vary  input  para- 
meters make  them  well  suited  to  laying  out  the  expected  acoustic 
ranges  for  any  area  in  a  style  best  suited  to  the  situation. 
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Greater  detail  and  lesser  roundoff  error  can  be  obtained  by  a 
proper  choice  of  display.   At  present  this  option  is  not  available 
to  ICAPS.   Also,  this  study  showed  that  for  short  time  periods, 
the  superiority  of  any  climatology  versus  another  is  unfounded. 
Yet  for  any  longer  period  of  consideration,  the  much  more  sophisti- 
cated climatology  of  FNWC  as  opposed  to  the  interim  climatology  of 
NAVOCEANO  should  present  more  representative  information. 
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V.   RECOMMENDATIONS 

A.  STANDARD  DEVIATION  FIELDS 

In  order  to  supply  the  ICAPS  operator  with  a  scheme  upon  which 
to  base  reliability  of  his  forecasts,  it  is  suggested  that  NAVOCEANO 
compile  Lewit's  (1971)  tables  in  such  a  way  that  the  expected 
standard  deviation  from  the  mean  temperature  field  for  each  acousti- 
cal area  can  be  made  operationally  useful.   A  color  code  for  large, 
moderate  or  small  standard  deviation  fields  is  a  simple  but  useful 
guide.   Such  a  presentation  will  afford  the  operator  an  insight  into 
how  much  variability  can  be  expected  in  his  area.   Also,  he  could 
plan  his  data  collection  based  on  such  a  system,  e.g.  more  BT's 
could  be  taken  in  areas  of  large  standard  deviation  and  fewer  in 
areas  of  small  standard  deviation. 

B.  TRAINING  OF  ICAPS  OPERATORS 

NAVOCEANO  and  the  ASW  schools  should  conduct  a  short  course  on 
the  proper  operation  of  the  ICAPS  system.   Ships  equipped  with  ICAPS 
should  be  assigned  trained  oceanographers  with  a  sense  of  apprecia- 
tion for  real-time  problems  and  the  art  of  good  data  collection  and 
forecasting  methods. 

C.  ALTERATIONS  TO  ICAPS  SYSTEM 

BT  input  should  be  recorded  on  a  cassette  tape  for  easy  handling 
and  recovery.   Presently  a  profile  must  be  hand  digitized  and  entered 
into  the  computer.   Also  the  program  should  be  altered  to  allow  for 
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the  input  of  more  than  a  single  set  of  parameters  at  a  time.   The 
program  could  solve  and  display  each  set  separately,  which  would 
save  considerable  time. 

D.   FNWC  PROGRAM  ALTERATIONS 

Presently  FNWC  calculates  a  variable  constant  which  is  an  indi- 
cation of  the  rapidity  with  which  the  mean  temperature  profile  may 
vary  from  its  initial  state  (Clark,  personal  communication) .   This 
constant  is  based  upon  the  data  density  (number  of  observations 
used  to  form  the  statistics) ,  the  sonic  layer  depth  and  the  magnitude 
of  the  sound  velocity  gradient.   The  number  is  then  simply  sent  to 
the  user. 

Two  basic  problems  are  noticed  in  this  approach.   Parameters 
such  as  expected  standard  deviation  of  the  temperature  profile,  and 
actual  roughness  and  slope  of  the  bottom  spatial  and  temporal  varia- 
tions in  received  BT  profiles  are  not  included  in  the  program.   More 
importantly,  the  user  does  not  know  which  parameter  exerted  the  major 
influence  in  the  determination  of  the  constant.   For  the  first  case, 
the  program  should  be  rewritten  to  include  the  three  parameters 
listed  since  they  exert  a  large  effect  on  the  variability  of  the 
acoustic  profile.   Secondly,  for  the  latter  problem,  the  user  should 
be  informed  as  to  which  parameter  had  the  greatest  influence  on  the 
forecast  in  order  to  enable  him  to  better  deal  with  the  problem.   A 
simple  code  could  be  established  showing  the  degree  of  relative  in- 
fluence each  had  on  the  acoustic  profile.   In  this  way  the  usefulness 
of  a  single  BT  is  enhanced  in  that  a  user  would  be  able  to  judge  as 
to  the  variability  of  the  water  column,  spatial  variation,  and  other 
parameters  separately. 
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Therefore,  it  is  the  opinion  of  this  investigator  that  the 
program  be  rewritten  to  include  additional  parameters  such  as 
those  discussed  above.   Additionally,  some  method  to  inform  the 
user  of  the  most  influential  parameter  should  be  developed. 
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APPENDIX  A 
FNWC  SEA  SURFACE  TEMPERATURE  ANALYSIS 

The  basis  for  the  Fleet  Numerical  Weather  Central's  (FNWC)  sea 

surface  temperature  analysis  scheme  is  the  Fields  by  Information 

Blending  (FIB)  program.   Figure  9  shows  a  flow  diagram  for  the  FIB 

program.   The  analysis  scheme  requires  a  first  guess  temperature 

field,  T  ,  and  a  base  field  temperature,  T  .   Figure  10  shows  the 
G  B 

area  module  for  the  program.   Gradients  are  given  as   u  and   v 

with  associated  weights  B*  and  C*  (Holl  et  al,  1971). 

The  basis  for  T  is  the  preceding  analysis,  T*   ,  where    t 
G  -T 

is  the  analysis  period,  adjusted  slightly  toward  climatology.   The 
climatology,  T  ,  is  computed  as  a  weighted  mean  for  each  analysis. 
The  weight  assigned  to  the  current  month  is  one-half, and  the  preced- 
ing and  following  months  have  a  combined  weight  of  one-half,  pro- 
rated by  the  hour  of  the  month.   Thus  the  equation  for  the  first 
guess  tenvperature  is  as  follows: 

T   =  (1-F) «T*  +  F'T 
G  -X  c 

where  F  is  proportional  weight  assigned  to  climatology  (Holl  et  al, 

1971). 

In  order  to  nullify  the  small  scale  anomalies,  often  caused  by 

sparse  data,  a  climatic  smoothing  is  accomplished: 

T  =  (T  -  T  )         +  T  . 
B      G     C  SMOOTHED     C 

The  first  guess  difference  field, J  ,  is  given  by 

*7   =  T   -  T  .  (Holl  et  al,  1971) 
v/  O     G     B 


67 


Figure  9.   Flow  diagram  for  SST  analysis 
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Figure   10.      Grid   diagram  for   SST  analysis 
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Observations  (subscript  n)  are  transformed  into  difference 
values: 

H       =  T   -  T  .  (Holl  et  al,  1971) 
<J  n    n    B 

Data  preparation  is  a  screening  process  designed  to:   (1)  elimi- 
nate all  reports  not  containing  an  SST  value;  (2)  eliminate  duplicate 
reports;  (3)  reject  reports  outside  analysis  region;  (4)  assign 
report  reliability,  A  ,  (a  function  of  the  contributing  variances) ; 
(5)  perform  an  initial  gross  error  check  by  analysis  the  magnitude- 
of  J     ;  and  (6)  reduce  reports  weights  for  quasi-duplication,  that 
is,  for  non-identical  reports  from  the  same  reporting  platform 
(Holl  et  al,  1971) . 

During  the  assembly  phase,  each  report  is  assembled  at  the  nearest 
grid  point.   The  report  difference  value,  ^/   ,  and  its  weight,  A  , 
are  added  to  the  accrued  information  at  the  grid  point: 

\  =  \  +  Vi 

Vn   =  Vn  +  Vl  Vl)/AN- 


N  refers  to  the  stage  of  the  assembly  and  no  to  the  single  report 
(Holl  et  al,  1971) . 

The  blending  equation  reduces  to 

^£,m=  S£,m[A£,m,N  £,m,N  +  {B£,m%  m+1  +Cl,m  l+ltm+Bl,m-l   lmm-1 


+co  -.   Vo    1       }'• 
£-l,m  £-l,m 


where  S  is  defined  by 


* 


S0   =  A0   +  B.    +  C0   +B„    ,  +  C„  . 
X-,n   £,n   £,n    x,,n   x,,m-l    £-l,m 


A*  is  solved  for  using  the  following  equation: 

*  _i    *   _1 

An    =  A.   /A.   /A. 

Ji,m    Z,m     £,m  £,m 
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If  this  is  the  first  pass  through  the  Assembly,  J*  Blending, 
and  A*  Solution  phases,  a  gross  error  check  is  now  made  on  each 
SST  observation.   All  reports  are  re-evaluated  considering  the  new 
J*  field.   Once  again  the  accepted  reports  are  assembled  at  the 
nearest  grid  point  and  the  J*  Blending  and  A*  Solution  phases  are 
completed.   This  time  the  output  is  combined  with  the  Base  Field 
Temperature  (T  )  to  provide  the  analyzed  sea  surface  temperature 
field  (Holl  et  al,  1971). 
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APPENDIX  B 
FNWC  MIXED  LAYER  DEPTH  ANALYSIS 

The  FNWC  mixed  layer  depth  (MLD)  analysis  scheme  is  described 
by  Hesse  (1973) .   The  MLD  is  defined  as  the  lower  boundary  of  the 
turbulent,  mixed  surface  layer  or  the  upper  boundary  of  the  thermo- 
cline.   The  potential  mixed  layer  depth  (POTMLD)  is  defined  as  the 
depth  to  which  the  turbulent  mixed  layer  depth  extends.   It  repre- 
sents the  mean  MLD  about  which  short-period  fluctuations  can  be 
expected. 

The  model  has  three  steps: 

1.  Adjusting  the  first  guess — the  12-hour-old  MLD  analysis — 
climatology  to  form  Guess  1. 

2.  Computing  the  effects  of  wave  action  on  the  MLD,  and 
adjusting  Guess  1  accordingly. 

3.  Adjusting  the  guess  field  (Guess  1)  to  available  BT  reports 
taken  in  last  72  hours. 

In  the  actual  program  computations  the  analysis  proceeds  as 
follows: 

1.   Guess   =  12-hour-old  POTMLD 

Guess  =  Guess   +  1/16  (CLIMAT  -  Guess  ) ,  where 

if  day  =  14,  CLIMAT  =  (0.5  -  day/30)  +  (day/30  +  0.5), 

if  day  =  15,  CLIMAT  =  CLIMAT  for  current  month, 

if  16  =  day  =  30,  CLIMAT  =  (1.5  -  day/30) (CLIMAT  current 
month  +  (day/30  -  0.5)  (CLIMAT  next  month), 

if  day  =  31,  CLIMAT  =  same  as  for  30. 
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2.  Compute  MLD: 

MLD       =  CH(11  -  0.1 (SST-  T,_n) ) ,  if  SST  =  Tc_- 
waves  600  600 

=  CH  •  11  if  SST  =  T   _, 

60U 

where  CH  =  combined  wave  height, 

SST  =  sea  surface  temperature, 

T  or,  =  temperature  at  600  feet. 
600 

3.  Smooth  the  guess  field  with  a  special  smoother. 

4.  Force  the  guess  field  to  be  at  least  as  deep  as  the  MLD 

waves  field: 

Guess   =  Guess, ,  if  Guess,  ~      MLD 

2        1  1       waves 

=  MLD      ,  if  Guess ,<  MLD. 
waves  1 

5„   Adjust  Guess   field  to  BT  reports. 

6.   Limit  final  analysis  field  by  forcing  POTMLD  =  900  feet. 
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Figure  HA.   Basic  flow  chart  for  sub-surface  temperature  analysis 
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Figure  11B.   Basis  flow  chart  for  sub-surface  temperature  analysis 

(continued) 


Read  in  GUESS  3  field  for  lev.L 


Read  in  DATA  LIST  for  level  L 


Perform  FIRST  Pass  tHrough  ANALYSIS  routine:  Form 
ANAL  1  by  adjusting  GUESS  3   to  reports  which 
passed  GEC,  using  L.   Carstensen's  analysis  tech- 
nique (XR);  reevaluate  data  by  cross-checking 
neighboring  reports 


s.'. 


Apply  VERTICAL  CONTROLS:  (1)  Apply  a  light  smoother 
to  the  vertical  gradient  between  ANAL  1   and  the  fi 
nal  analysis  at  the  level  above.  (2)  Require  T(L) 
T  (level  above)  except  at  selected  points.  ANAL  2 


_s' 


Perform  SECOND  PASS  through  ANALYSIS  routine: 

Form  ANAL  3  by  adjusting  ANAL  2  to  data,  using  XR 
with  lateral  check 


Apply  VERTICAL  CONTROL:  require  T(L)  T( level 
above)  except  at  selected  grid  points. 
Result  is  ANAL  4. 


Restrict  T(L)  to  range  (-)2  C  to  (+)35  C. 
This  provides  final  analysis  for  level  L. 


L  =  next  level 


75 


APPENDIX  C 
FNWC  SUB- SURFACE  TEMPERATURE  ANALYSIS 

The  FNWC  sub-surface  temperature  analysis  is  run  at  the  100, 
200,  300,  400,  600,  800  and  1200  foot  levels.   The  12-hour  old 
analysis  is  used  as  a  starting  point  for  each  level.   This  is 
adjusted  1/16  of  the  way  toward  climatology .   After  this  adjust- 
ment, the  figures  are  further  adjusted  to  the  MLD  and  SST  analy- 
sis to  produce  the  guess  field.   The  MLD  and  SST  analysis  serve 
as  strong  anchor  points  for  the  subsequent  sub-surface  analysis. 
(Hesse,  1973) 

BT  reports  for  the  preceding  72-hour  period  are  included  in 
the.  analysis,  and  they  are  weighted  by  age  on  a  1  to  4  scale. 
At  each  level  the  report  is  subjected  to  a  gross  error  check  (GEC) 
by  requiring  that  the  report  fall  within  limiting  values  of  the 
guess  field.   A  similar  GEC  is  run  on  reported  gradients  by  com- 
parison with  the  guess  field  gradients.  (Hesse,  1973) 

The  guess  field  is  adjusted  to  all  reports  which  passed  the 
GEC.   The  vertical  gradient  between  the  first  analysis  and  the 
final  analysis  from  the  level  above  is  then  slightly  smoothed. 
Next,  consistency  is  checked  by  allowing  a  positive  temperature 
gradient  to  exist  only  at  selected  grid  points.   Finally,  the  tem- 
perature range  at  all  levels  is  restricted  to  the  range  of  -2 
degrees  to  3  5  degrees  centigrade.   Figure  10  shows  the  basic  flow 
chart  for  the  sub-surface  analysis.   (Hesse,  1973) 
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